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Fermionic natural occupation numbers do not only obey Pauli’s exclusion principle, but are even further re¬ 
stricted by so-called generalized Pauli constraints. Such restrictions are particularly relevant whenever they are 
saturated by given natural occupation numbers A =s (A*). For few-site Hubbard models we explore the occur¬ 
rence of this pinning effect. By varying the on-site interaction U for the fermions we find sharp transitions from 
pinning of A to the boundary of the allowed region to nonpinning. We analyze the origin of this phenomenon 
which turns out be either a crossing of natural occupation numbers A i(U), A*+i (U) or a crossing of iV-particle 
energies. Furthermore, we emphasize the relevance of symmetries for the occurrence of pinning. Based on 
recent progress in the field of ultracold atoms our findings suggest an experimental set-up for the realization of 
the pinning effect. 
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I. INTRODUCTION AND CONCEPTS 

The Pauli exclusion principle restricts fermionic occupation 
numbers rii , 

0<n*<l. (1) 

From this 1-particle viewpoint the strong influence of Pauli’s 
exclusion principle for fermionic quantum systems is not sur¬ 
prising. Since it prevents fermions in lattice systems from 
hopping to a neighboring site which is already occupied by 
another fermion in the same spin state it can restrict the mo¬ 
bility of fermions significantly. 

It has been suspected since the 1970s Qi but was shown 
only recently 113-5] that the fermionic exchange statistics 
does not only imply Pauli’s exclusion principle but leads 
to even stronger restrictions on fermionic natural occupa¬ 
tion numbers. To be more specific, we consider N identi¬ 
cal fermions described by pure antisymmetric quantum states 
|\P) G A n [H^]. Here, T~L^ is the underyling d-dimensional 
1-particle Hilbert space. To each |\P) we can assign its natural 
occupation numbers (NON) A i, the eigenvalues of the corre¬ 
sponding 1-particle reduced density operator 

d 

p 1 =NTr ff - 1 mm = '52 X k\ k )(k\- ( 2 ) 

k= 1 

pi is obtained by tracing out N — 1 fermions and its eigen¬ 
states, | k), are called natural orbitals. Note, that the implicit 
dependence of A& and | k) on pi and |\P), respectively, is sup¬ 
pressed. Pauli’s exclusion principle © can be reformulated 
as 

0 < A, < 1. (3) 

The so-called generalized Pauli constraints (GPC) then take 
the form of linear inequalites 0-0] 

Df' d \ X) s «<“> + £««>A, >0, (4) 
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(i) 

with n - G Z, j = 1,..., < oo, which give rise 

to a poly tope VN,d C of allowed vectors A = (A i)f = i 
of decreasingly-ordered NON. The family of GPC implies 
Pauli’s exclusion principle and V^,d is in particular a proper 
subset of the ‘Pauli hypercube’ [0, l] d described by (O. 

A direct significance of GPC was proposed by Klyachko 
16, 7] for ground states: For some systems the minimization 
process of the energy F7[\I/jv] = could get stuck 

on the boundary of the poly tope Vna since any further mini¬ 
mization would violate some GPC. In that case, ground state 
properties would not only rely on the form of the Hamiltonian 
but also significantly on the structure of GPC. Such pinning of 
A to the polytope boundary leads to further remarkable conse¬ 
quences. On the one hand, it potentially restricts the dynamics 
from the 1-particle viewpoint since A can never leave the poly¬ 
tope. On the other hand, it implies that the corresponding N- 
fermion quantum state has a significantly simplified structure 
and also reduced entanglement Bn®. 

GPC may have a more indirect physical relevance, as well. 
For instance, they may lead to improvements in reduced den¬ 
sity matrix functional theories as explored in 111,112], can be 
used for describing the openness of fermionic quantum sys¬ 
tems ill and allow to quantify the influence of the fermionic 
exchange statistics beyond that of Pauli’s exclusion principle 
01 . 

Based on these ideas, a fundamental question arises: Do 
there exist fermionic systems exhibiting the pinning-effect? 
Although Klyachko suggested a mechanism for pinning it 
would be rather surprising if NON of interacting fermions did 
exactly saturate some of those 1-particle constraints. For in¬ 
stance, NON for interacting fermions do never saturate Pauli’s 
exclusion principle, since none of them is ever exactly iden¬ 
tical to 1 or 0. Further evidence that pinning is unlikely to 
occur is provided by reduced density matrix functional the¬ 
ory. It seeks a distinguished functional F[pi\ whose mini¬ 
mization leads to the energy and the 1-particle reduced density 
operator pi of the ground state 03. Since ^[pi] is nonlinear 
(see e.g. 116]), there is no reason why the minimum of ^[pi] 
should be attained on the polytope boundary. 

Given all this evidence against pinning, it is quite surprising 
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that its occurrence was claimed in jhjfor the beryllium atom 
based on numerical data provided in J17l]. However, by taking 
additional digits into account (in [17|] only the first six were 
shown), this can be disproved Hill . Moreover, in [19] first 
analytic evidence was found that NON for ground states of 
interacting fermions do lie indeed close to (but not exactly on) 
the polytope boundary. The occurrence of this quasipinning 
[ 8, fg] was also found for several small atoms and molecules 
Ii20l - l22ll . provided that the numerical approximations were not 
too restrictive (23]. All these observations suggest that the 
occurrence of pinning for non-approximated ground states is 
highly non-generic, or even impossible. 

In this paper, we provide in form of the few-site Hubbard 
model a first example for a system exhibiting pinning. Be¬ 
sides its importance for solid state physics the Hubbard model 
has gained much relevance in recent years on the microscopic 
scale, as well. This is due to progress in the field of ultra¬ 
cold fermionic gases which allows to study the crossover from 
few-fermion to many-fermion physics (see e.g. 0, 0). In 
that context, the Hubbard model was realized experimentally 
very recently [126]. The corresponding Hamiltonian (in second 
quantization) reads 


r—1 r —1 

H = ~t ( C Ula C ia + h ' C ') + U J2 "4 ( 5 ) 

2=0 a z=0 

where c\ a and Ci a are the fermionic creation and annihilation 
operators for a spin— fermion on the i- th lattice site with spin 

a , l with respect to the z -axis and hi a = c\ a Ci a . 

The paper is organized as follows. In Sec. [TH we recall the 
symmetries of © and introduce a symmetry-adapted basis of 
TV-fermion quantum states. The main results are presented in 
Sec. mu There, we analytically determine the eigenstates of 
© for three fermions on three lattice sites and explore possi¬ 
ble pinning. This is extended by an exact numerical approach 
for the next larger settings involving more fermions and/or 
more lattice sites in Sec. [TV] In Sec. [Vj we explain the role 
of symmetries for the occurrence of pinning. By building on 
very recent progress in the field of ultracold gases we propose 
in Sec. |VI|two ideas for experimental realization of pinning. 
Sec. lYIIl provides a short summary and a conclusion. 


II. SYMMETRIES 


In this section we present the most elementary symmetries 
of the Hubbard model © and introduce a basis of symmetry- 
adapted quantum states. 

Since the symmetries of the Hubbard model are well- 
known (see e.g. Ref. they are just listed in order to 

keep our paper self-contained. First of all, since the Hamil¬ 
tonian H commutes with the total fermion number operator 
N = Y^=o Sa ^ 2 , [H, TV] = 0, the total particle num¬ 
ber N is conserved and we restrict to fixed fermion numbers 
N. By employing first quantization in the following H is 
then restricted to the corresponding TV-fermion Hilbert space 
Hn = of antisymmetric quantum states, where the 


2r-dimensional 1-particle Hilbert space has the substructure 
Hi = H^ (8) H{ s \ Here, H^ = C r describes the spatial (/) 
degrees of freedom and h[ s ^ = C 2 is the spin Hilbert space 
of a single fermion. 

We denote the total spin vector operator by S, its z- 
component by S z and let T be the translation operator which 
translates each of the N fermions from its lattice sites i to the 
next site i - hi. Then, the symmetries of H are described by 
the following relations (see e.g. 0) 

[hJ 2 ] = 0 , [H,S z ]=0, [H,T] = 0 . ( 6 ) 

The operators S 2 ,S Z and T also commute with each other. 
Consequently, the total spin quantum number S, the magnetic 
spin quantum number M and the total Bloch number (wave 
number) K are good quantum numbers ll28ll . The Hamilto¬ 
nian © is block diagonal w.r.t. those symmetries and the total 
Hilbert space splits according to 

T~Ln = A n [H^ <S>h[ s ^} = (J) (J) (J) Hs,m,k • (7) 

S=S- M= — S K =0 

By setting h = 1 the maximal total spin is given by S+ = ^ 
and the minimal total spin S- by 0 for N even and \ for N 
odd. 

Since we will focus below on the 1-fermion picture it is 
worth noticing that the 1-particle reduced density operator 
pi inherits symmetries from the corresponding TV-fermion 
quantum states |\h) 0. Whenever |\h) has a symmetry, 
F® |\h) = e^|\k), where F\ is a 1-fermion operator act¬ 
ing on the corresponding pi inherits the same symme¬ 
try, [pi, Fi] = 0 and is therefore block-diagonal w.r.t. the Fi~ 
eigenspaces. 

For the Hubbard model with periodic boundary conditions, 
the lattice translation T = T® defines such a symme¬ 
try, where Ti is the 1-fermion translation operator. Another 
symmetry is generated by the ^-component of the fermion 
spin. Consequently, the natural orbitals arising from any 
symmetry-adapted state |\h), |\h) G 1~Ls,m,k, are given by 
| kcr) = |fc) (8) |cr). Here | k) G \ k = 0,1,..., r — 1 is the 
corresponding Bloch state, Ti\k) = e^ lk \k), and |cr) G 
, cr ="[, is the corresponding spin state. Moreover, since pi 
is diagonal w.r.t. the NO, we can easily calculate the NON. 
They are given by the diagonal elements (ka\pi\ka). 

The basis {| kcr)} for Hi induces a basis for the TV-fermion 
Hilbert space Hat, the TV-fermion Slater determinants 


\ki(Ji ,..., k N a N ) = A n [\ki(Ji) ( 8 )... ( 8 ) \k N a N }} , (8) 

where An is the A/-particle antisymmetrizing operator and 
the ordering of the 1-particle basis states is given by |0 t), 
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III. THREE FERMIONS ON THREE SITES 

In this section we study the Hubbard model ([5]) for three 
fermions on three lattice sites with periodic boundary condi¬ 
tions 0 . We analytically diagonalize the Hamiltonian and 
explore whether its eigenstates exhibit the pinning or quasip- 
inning effect. 


A. Analytic diagonalization of the Hamiltonian 

To diagonalize (0 on the corresponding ( 3 ) = 

20-dimensional 3-fermion Hilbert space we first block- 
diagonalize H w.r.t. to its symmetries which were mentioned 
in Sec. [Ill The total spin quantum number S can take the 
two values S = §, We first split the total Hilbert space 
w.r.t. those two values for S and then continue with the mag¬ 
netic quantum number M : 

• S = §. The only state with maximal magnetic quan¬ 
tum number M — | is the symmetry-adapted 

state |0 t?l t)- It has the Bloch number 
K = 0 + 1 + 2 (mod 3) = 0. By successively 
applying the lowering ladder operator to that state 
we find the remaining three symmetry-adapted 
states of the corresponding quadruplet. They 
read [\0i, 4, 2 f> + | 0 f, U, 2 f> + | 0 f, It, 2 t)], 

75 [|0t, 4, 2t> + |0t, 4 , 2t> + |0t, 4 , 2t>] and 

|0t, It, 2t). As an elementary exercise one verifies that 
© restricted to that 4-dimensional subspace is given 

by-ff|s=| = 0 . 

• S = The corresponding subspace is 16-dimensional. 
Since the Hamiltonian is also invariant under simultane¬ 
ous flips of all spins w.r.t. the z- axis, we restrict without 
loss of generality to M > 0, i.e. M — 

• M = In this 8 -dimensional subspace we 
find two states with K = 0, [|0 11, 2 t) 

+ jO t, U, 21) + |01,11, 21)], and 

the Hamiltonian restricted to that 2-dimensional 
subspace, H 1 1 0 , is given by u I 2 . The remain¬ 
ing subspaces are H 1 1 x and Hi i 2 - However, 

due to the additional invariance of H under inver¬ 
sion in the reciprocal lattice (recall that K = 2 
equals K = —1), we can restrict to the case 
K = 1. This subspace is spanned by the three 
states 101 , 01 , 11 ), 111 , H, 21 > and 101 , 21 , 21 ) 
and H represented w.r.t. to those three states takes 
the form 



The diagonalization of this 3 x 3-matrix leads to 
a cubic equation which is presented and solved in 
Appendix [A] 


• M = — \. Results are identical to the case M = \ 
up to a simultaneous flip of all three spins w.r.t. the 
z-axis. 

This completes the diagonalization of the Hubbard Hamilto¬ 
nian ( 0 . 

The energy spectrum of H is shown in Fig. 0 The three 



u 


FIG. 1. Energy spectrum for the 3-site Hubbard model for three spin- 
\ fermions. Each eigenvalue has multiplicity four. The dashed black 
line (horizontal line) corresponds to S — §. The case S — \ yields 
the black straight dotted line (K = 0) and the three red solid lines 
correspond to the solution for the nontrivial subspace problem ([9]) for 
K = ± 1. 

solid lines E\ (u ), E 2 (u) and E 3 (u) are those arising from the 
eigenvalue problem for ©, 

H\hi ll |^(«)) = s # («)|^(u)), j = 1,2,3. (10) 

2 ’ 2 ’ 

According to the symmetries each Ej (u) has multiplicity four 
(quantum numbers S = M = ±| and K = ±1). For the 
four eigenstates with S = § we find the eigenvalue E{u) = 0 
(dashed line) and for those with S = \,M = and K = 0, 
E(u) — u (dotted line). 


B. Pinning analysis 

In this section we investigate whether the energy eigen¬ 
states of the Hubbard Hamiltonian are showing quasipinning 
or pinning. Due to their particular relevance for experiments 
we first consider symmetry-adapted states. At the end of this 
section we also study coherent superposition of degenerate 
eigenstates. 

The corresponding GPC (cf. Eq. (0) for the setting of N = 
3 fermions and a 6 -dimensional 1-particle Hilbert space are 
given by d] 

Ai + A6 = A2 + A5 = A3 + A4 = 1 
2)( 3 ’ 6 )(a) = A5 + Ae — A4 + 0 . (11) 

The first three GPC lead to strong structural implications for 
the corresponding 3-fermion quantum state (see e.g. 0 ). 
Since they are always saturated it only makes sense to explore 
a possible saturation of the fourth constraint, which takes the 
form of a proper inequality. 

First, we study the four eigenstates with S = |. Those 
two with a minimal/maximal M are single Slater determi¬ 
nants and consequently their NON equal the Hartree-Fock 
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point A hf = (1,1,1,0,0,0). Since this vector even satu¬ 
rates the weaker Pauli exclusion principle constraint dTJ it is 
also trivially pinned to the polytope boundary. The other two 
eigenstates, those with M = ±^, have NON |(2, 2, 2,1,1,1) 
which is not pinned. This example also demonstrates that 
NON of degenerate energy eigenstates can be quite different. 

Now, we continue with the S = ^-eigenstates. Those four 
with K — 0 also do not depend on u and their NON are given 
again by | (2, 2, 2,1,1,1). More interesting are the remaining 
three eigenstates (with multiplicity four) arising from (0 and 
m . since they do depend on u. These states have the general 
form 

Mu)) = lot, 0|, It) + Pin) I It, 4, 2t> 

+ 7 (ix)|0t,2t,2;). (12) 

The three coefficients a(u),/3(u) and 7 (u) are calculated in 
Appendix lAl and their absolute squares are presented in Fig. [2] 
for the three non-trivial eigenstates following from (ITOt . 


Furthermore, it is an elementary exercise [§] to show that 
the normalization of p\ and p\ together with the first three 
GPC (fTTI) implies that the largest two NON, Ai and A 2 , belong 
to the block p[, the smallest two NON, A 5 and A 6 , belong to 
the block p\ and the smallest eigenvalue in of p\ and the 
largest eigenvalue of p\ fulfill 


Eq. (fl5l) simplifies the pinning analysis and we can distin¬ 
guish two cases 

• n rnax — n min : The fourth-largest NON A 4 is therefore 
given by in and for the second line of Eq. (fill) we 
find 

Z)( 3,6 )( A) = A5 + A6 — A4 

(1 71max ) ^min ^ ’ (^) 



where we used that A 5 and A 6 both belong to p\, 
Tr [p\] = A 5 + A 6 + nj^ ax = 1 and Eq. (Q3]). 

• n rnax < n lnin • I n a similar way we find for this case 


D ^ 6 \ X ) = ( l - nt ax )- n , 


max 

— n f — rr f 0 

mm max ^ w • 


(17) 


FIG. 2. The absolute squares of the coefficients a(u) (red), /3(u) 
(blue) and 7 (u) (green) for the eigenstates following from (fTol . 
Those for the ground state |^i (u)) (left) are identical to those for the 
highest energy state |^ 3 (^)) (not shown) up to a reflection w.r.t. the 
vertical axis. On the right the corresponding result for the excited 
state 1 4*2 (^ 0 ) shown. 

The corresponding 1-particle density operator represented 
w.r.t. to the 1 -particle states |fccr) is diagonal as explained in 
Sec.UH In particular, this means that the two different spin- 
blocks are decoupled [3l]], i.e. 

Pi = Pi® Pi ■ (13) 

Since pi is normalized to TV = 3 and needs to reproduce the 
5^-spin expectation value M = \ the blocks have the specific 
normalizations Tr[p|] = 2 and Tr[pf ] = 1. 

For the eigenvalues k = 0,1, 2 of p\^ corresponding 
to the natural orbitals |fc, t / i) we find 

n l = H u )\ 2 + It ( m )| 2 A = \ a ( u )\ 2 
n \ = \ a ( u )\ 2 + \/ 3 ( u )\ 2 , n \ = l /3(«)| 2 
A = \P{u)\ 2 + |7(w )| 2 A = |7 («)| 2 • (14) 

The structure of this spectrum - three eigenvalues are given 
by the weights |ck| 2 , \/3\ 2 and | 7| 2 and the other three by 
sums of two of them - is a consequence of the symmetries, 
only. The concrete form (beyond the symmetries) of the Hub¬ 
bard Hamiltonian affects only the weights \a(u)\ 2 , \/3(u)\ 2 
and | 7 (^)| 2 . 



FIG. 3. u-dependence of the NON A i (u) for the ground state 147 (u)} 
(left panel) and for the excited eigenstate |^ 2 (^)) (right panel). The 
NON for the ground state are identical to those for the highest energy 
state 1 4 /3 (u-)) (not shown) up to a reflection w.r.t. the vertical axis. 
The dashed lines are the NON corresponding to spin f and the solid 
lines correspond to f. 

This means, depending on the ratio of and n^ n , we 
either have pinning or nonpinning. Whether the eigenstates 
(IT2b are pinned or not is illustrated in Fig. [3] It shows the 
six NON for the relevant on-site interaction regime. The three 
eigenvalues of the block p\ are shown as solid lines and those 
for p\ as dashed lines. Pinning is given for those ^-regimes 
where > n^ n , i.e. whenever the highest dashed line is 
above the lowest solid line. E.g. for the ground state |fl/i («)) 
presented on the left side of Fig.[3]this is the case for v < «o ~ 
12.86. At that ‘critical’ point, uq, there is a transition from 
pinning to nonpinning according to the crossing of (u) 
and n^^u) For the excited state ^(rO) presented on the 
right side we can expect an even more complex behavior, since 
the relevant NON curves do cross several times. 
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FIG. 4. The eigenstates dTot of the Hubbard model for three fermions 
on three sites show a pinning-nonpinning transition. The distance 
L^ 3,6 ) (u) of the NON A (u) to the polytope boundary for the ground 
state |T^i(u)) (top) is identical to that for the highest exited state 
|T^ 3 (u)) (not shown) up to a reflection u —u. The excited state 
|T^ 2 (u)) (bottom) shows a reentrance behavior of pinning in form of 
four pinning-nonpinning transition. 


Indeed, as shown on the upper side of Fig. @1 D^ 3 , 6 \u) 
for the ground state undergoes a pinning-nonpinning transi¬ 
tion. D( 3 , 6 \u) for the excited state |T r 2 (rt)) (lower panel of 
Fig.0]) shows indeed that an expected re entrance phenomenon 
is present, a sequence of pinning-nonpinning transitions. 

Due to the distinguished role of the ground state we de¬ 
termine its leading behavior on the right side of the ‘critical’ 
point uq. From the concrete form of the ground state |\Fi M> 
(see Appendix lAl) one obtains 

D^ 3,6 \u) ~ 0.026 (u — uo) , foru>uo. (18) 

So far we have restricted ourselves to symmetry-adapted 
eigenstates. However, since experimental realizations of the 
Hubbard model typically prepare states with fixed spin quan¬ 
tum numbers (see e.g. (26]) the most general ground state 
takes the form 

l*€,c(«)> = £l*i(«)) + CXqI*i(«)> , (19) 

but has still well-defined spin quantum numbers 5 = |, 
M = ±|. Here |£| 2 + \(\ 2 = 1, Tq is the 3-fermion in¬ 
version operator for the reciprocal (Q) lattice and |4>i (u)) has 
the form (IT2l) where the coefficients follow from the eigen¬ 
value problem © and (ITOb . We set ( = \C,\e lip and without 
loss of generality we choose £ = |£| = yj 1 — |C| 2 . Due to 
a ^-symmetry of the NON (see Appendix [B]) we present the 
pinning results in Fig. [5] for the two extremal cases p = 0, 



FIG. 5. Saturation D(u) for the ground state |^(u)) = 

VI — r 2 \ i (u)) +r e lip lQ \ 1 (u) ) of the Hubbard model with three 
sites and three fermions for the two extremal cases p = 0 (left) and 
P = f (right). 


There we can see that pinning is uniform in u whenever 
ICI 2 = |£| 2 = r 2 = For all other superpositions it vanishes 
for some critical value uq(£). Moreover, since this critical 
value uo(C) is minimal for either ( = 0 or \(\ = 1 , we con¬ 
clude that superposing |\Fi(rO) an d ^q|^i(^)) extends the 
ti-regime of pinning from (—oo, 12 . 86 ] to some larger inter¬ 
val (—oo, ixo(C)]* u o(C) ^ 12.86. Furthermore, as already 
suggested by Fig. [5] it turns out that the relative phase p = \ 
favors pinning best. 


IV. LARGER SETTINGS 


In this section we study the next larger systems. However, 
since complete families of GPC are not known yet for cases 
corresponding to more than N = 5 fermions and more than 
r = 5 lattice sites we restrict to N < 5 and r < 5. The 
GPC for all those cases can be found in & We diagonal¬ 
ize the corresponding Hamiltonians exactly by numerical 
methods. In the following we present as representative 
cases the two cases (TV, r) = (3,4), (5, 5). The other cases 
(3, 5), (4,4), (4, 5) are qualitatively similar. Note also, that 
due to the particle-hole symmetry (2711 . the NON for the cases 
(AT, r) with N > r follow from those for N < r. 


A. Three fermions on four sites 

For the Hubbard model with three spin— fermions on four 
lattice sites we explore possible pinning for several symmetry- 
adapted eigenstates. In Fig. [ 6 ] we present the results, i.e. the 
minimal distance D min (u ) of the vector A (u) of NON to 
the polytope boundary dTz^, for the 5-dimensional subspace 
with quantum numbers (5, M, K) = (|, |,1). The eigen¬ 
states for the other subspaces show similar behavior. D m i n (u ) 
for the five eigenstates are ordered from the front to the back 
according to increasing energy. The first and fifth as well as 
the second and fourth state behave identically up to a reflec¬ 
tion of u at u = 0 ill. D m i n for the eigenstate with minimal 
energy shows transitions between pinning (shown in black) 
and quasipinning, which is here defined as < 10 -3 . 

The second eigenstate in this subspace exhibits pinning for 
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negative u. For positive u there is a pinning-quasipinning- 
pinning transition followed by nonpinning for larger u. The 



third eigenstate is pinned in the interval [—3.75,3.75] and 
nonpinned in the remaining regime. 



u 


FIG. 6. For the Hubbard model with three fermions on four sites all five eigenstates with quantum numbers (S,M,K) = (|, 1) are 

explored w.r.t. possible (quasi)pinning. The curves Dmi n (u) showing the minimal distance of the NON A(u) to the polytope boundary are 
ordered according to increasing energy from the front to the back. Pinning is shown in black and the strength of quasipinning is visualized in 
form of different colors. 


It is instructive and also relevant for possible experimen¬ 
tal realizations of the pinning effect to study the global 
ground state separately. For the ^-interval (—18.6,18.6) 
the ground state belongs to the blocks (S,M,K) = 

(|, ±j>, ± 1 ) which were already discussed and presented in 
Fig. [ 6 ] At ‘critical’ points u ~ ±18.6 energy eigencurves 
cross and the ground state changes its quantum numbers to 
(S', M, K) = (±,±±,± 2 ) for u < -18.6 and (S,M,K) = 
(§, ±§, 0), (§, ±~, 0) for u > 18.6. The ground state pin¬ 
ning behavior is shown in Fig. [71 The behavior in the regime 
(—18.6,18.6) is already known from Fig. [ 6 ]and for \u\ > 18.6 
there is always pinning. The nonanalytic behavior of the dis¬ 
tance D min (u ) of the NON A (u) to the polytope boundary 
at ‘critical’ points has different origins. First, recall that the 
Hamiltonian ([5]) depends analytically on the on-site interac¬ 
tion u and the energy spectrum and eigenstates inherit this 
analyticity. However, the ordering of the eigenstates accord¬ 
ing to increasing energies may violate this analyticity. This 
happens as already mentioned above at the ‘critical’ point 
u = —18.6. At the other ‘critical’ point u = ±18.6 Dmin^u) 
is analytic since both quantum states whose energy curves are 
crossing show identical pinning behavior (both are pinned). 
Even for a given analytic 7V-fermion state |\f r (i/)) nonanalytic 
behavior of D min (u ) is possible. Although pi (recall ©) and 
therefore also its spectrum (NON) and natural orbitals inherit 
this analyticity the ordering of the NON A i(u) > A i+i(u) 
may violate it and could therefore lead to ‘critical’ points for 
Dmin{u). This happens at u ~ 2.1 (see inset of Fig. [7]). 
To verify this we study the distances of successive NON and 
present those three pairs in Fig. [ 8 ] which indicate crossings of 
NON. There is indeed a crossing at u ~ 2.1 since A 6 — A 7 
vanishes there. But why do the crossings at u = 4.4,12.1 
indicated in the first two graphics in Fig. [8] not lead to ‘crit¬ 
ical’ points for D m i n (u)l The reason is obvious. Crossing 



FIG. 7. For the ground state of the Hubbard model with three 
fermions on four sites possible (quasi)pinning is explored. For very 
negative on-site interactions u the ground states exhibits always pin¬ 
ning which vanishes at a ‘critical’ point u\ ~ —18.6 in form of 
a sharp pinning-nonpinning transition. Also quasipinning-pinning 
transitions can be found, namely at the ‘critical’ points U 2 = 0 and 
U 3 ~ 2.1 (see inset). 


of NON has only an influence on (u) if the correspond¬ 
ing two NON do enter D m i n differently (recall the form 0), 
i.e. with different coefficients K n . To verify this for the ground 
state studied in this section we first need to know the polytope 
facet which is closest to A (u), i.e. we need to find the most 
saturated GPC @. For u > 2.0 this is the constraint 

D{ A) = 2 — (Ai ± A2 ± A4 ± A7) ± 0. (20) 

Indeed, the NON A 6 , A 7 have different weights, = 0 and 
kj = —1. This is not true for the other pairs of NON, Ai, A2 
and A 5 ,A 6 since = ^2 = — 1 and ^5 = kg = 0. Con¬ 
sequently, the pinning behavior is also analytic around u ~ 
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FIG. 8. For the ground state of the Hubbard model with three 
fermions on four sites only the distances A i(u) — A*+i for i — 1, 5, 6 
of decreasingly ordered NON A i(u) have nonanalytic points. Such 
nonanalytic behavior for some ‘critical’ values for u can change the 
underlying pinning behavior Dmin ( u ) of the ground state as shown 
in Fig. [7] 


4.4,12.1. A third reason for nonanalytic behavior of D^^u) 
is given when the GPC which is most saturated changes, 
i.e. we have a crossing of two saturations D(u),D r (u). This 
happens at the point u & 2.0 (see sharp kink in the inset of 
Fig. 171) where the GPC (f20l) and D'(X) = 1 — Ai — A2 + A3 > 0 
interchange their role as most saturated GPC. 


B. Five fermions on five sites 

As a third case we study five fermions on five lattice sites. 
The ground state (with multiplicity four) belongs to the sub¬ 
spaces (S',M, K) = (^,± 7 p±l). Therefore, we explore 
possible (quasi)pinning for the eigenstates with (S, M.K) = 
(|, |, 1). The results are shown in Fig. [9j The eigenstates 
are ordered from the front to the back according to increasing 
energy. In particular, the first curve D m i n ( u) describes the ab¬ 
solute ground state. For several eigenstates pinning of NON A 
occurs only for the isolated point u = 0. In the neighborhood 
of u = 0 we find quasipinning which eventually changes to 
nonpinning for \u\ > 0.5. The eigenstates of the other sub¬ 
spaces with S = | show very similar behavior as those for 
K = 1 shown in Fig. [9] For S = §, f the corresponding sub¬ 
spaces are much smaller and the pinning behavior shows more 
transitions similar to the results for S = \ for three fermions 
on four sites presented in Fig.[ 6 l 

Finally, it should be also mentioned that linearly superpos¬ 
ing degenerate eigenstates belonging to subspaces with differ¬ 
ent quantum numbers leads again to enhancement of pinning 
as discussed at the end of Sec. lIIIBl 


V. ROLE OF SYMMETRIES AND STRUCTURAL 
IMPLICATIONS OF PINNING 

In this section we emphasize the strong relation between 
symmetries and pinning. One of the most remarkable conse¬ 
quences of pinning as an effect in the 1 -particle picture is that 
it allows to reconstruct the structure of the corresponding N- 
fermion quantum state |4/). In addition, |\F) is significantly 
simplified (see e.g. 1 6 , El-El])- To explain this, notice that for 
NON saturating all Pauli exclusion principle constraints, |\F) 



u 


FIG. 9. For the Hubbard model with five fermions on five sites all fif¬ 
teen eigenstates with quantum numbers (S, M, K ) = (|, |, 1) are 
explored w.r.t. possible (quasi)pinning. The curves D m in (u) show¬ 
ing the minimal distance of the NON X(u) to the polytope boundary 
are ordered according to increasing energy from the front to the back. 
Pinning is given (for each eigenstate) only for the single point u = 0 
(not visible in the graphic) and the strength of quasipinning is visu¬ 
alized in form of different colors. Quasipinning is found for several 
eigenstates in the neighborhood of u = 0. 


can be written as a single Slater determinant, 

A = A hf = (1,..., 1,0,..., 0) =^> |4>) = |ii,..., ijv) • 

N d-N 

( 21 ) 

This structural simplification for pinning to the Hartree-Fock 
point A hf, which is a vertex of the poly tope Vn^, general¬ 
izes to arbitrary points on the poly tope boundary: By expand¬ 
ing |\F) in Slater determinants built up from its own natural 
orbitals | i), 

W = ^ ‘ • • ^n) , ( 22 ) 

many of the expansion coefficients need to vanish 

rigorously in case of pinning. This selection rule of Slater 
determinant for pinning D( A) = 0 reads (see e.g. 110]) 

D\&\i) 7 ^ 0 => = 0. (23) 

Here | i) = |ii,..., ijv) and 

Dy = L>(ni, ...,n d ), (24) 

where fij is the particle number operator for the j-th natural 
orbital of |\F). Since D(-) is linear (recall Eq. d4j), D^ is a 
(hermitian) 1 -particle operator. 

The structural simplifications of 14/) by pinning lead to a re¬ 
duction of entanglement described, e.g., by the von Neumann 
entropy of the 1-particle reduced density operator p \. This fol¬ 
lows from the fact that the maximal entanglement corresponds 
to NON A = (2[,...,2jr), which never saturate any GPC. It 
is a future challenge to derive optimal upper bounds on the 
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FIG. 10. On the left side, the Pauli exclusion principle prevents the 
valence electron from decaying to lower lying energy shells. This 
kinematical restriction for A = (1,..., 1, 0,...) (E dV (red dot) 
generalizes to arbitrary ‘points’ A on the poly tope boundary dV 
(schematic illustration on the right side). 


entanglement entropy given pinning. Furthermore, since ex¬ 
citations depend qualitatively on the structure of the ground 
state, pinning may also influence the low temperature proper¬ 
ties. 

For systems of interacting fermions one does not expect 
any of the coefficients c^ of the superposition (l22l) to van¬ 
ish. Consequently, pinning seems unlikely or even impossible. 
However, as the analysis of the Hubbard model has shown, 1- 
particle symmetries, like the translational or spin symmetry, 
cause such a significant reduction of Slater determinants (see 
e.g. (fl2l) ) and eventually can cause pinning. Yet, our results 
also show that the presence of symmetries does not automati¬ 
cally imply pinning. 

Besides the main insight that symmetries of the Hamilto¬ 
nian favor pinning, a converse statement is also possible. Ac¬ 
cording to Eqs. (I22lh(l23l) and (l24k the occurrence of pinning 
reveals a 1-particle symmetry of the corresponding quantum 
state |\F). It is generated by the ‘observable’ D^ and reads 

e —|^f) = |^) (25) 

All consequences of pinning presented in this section hold 
approximately for quasipinning, as well 


VI. EXPERIMENTAL REALIZATION OF PINNING 


In this section we outline two conceptually different ideas 
for realizing and verifying pinning. Both of them are based on 
very recent progress in the field of ultracold fermionic gases 
simulating the few-site Hubbard model with full control over 
its quantum stateJH]. By focusing laser beams with a high- 
resolution objective, a potential landscape with a few wells 
is generated which can be loaded with 6 Li atoms. This sys¬ 
tem then simulates the few-site Hubbard model. By indepen¬ 
dently controlling the intensity and position of the laser beams 
with an acousto-optic deflector the hopping between the wells 
can be tuned[26]. Hence, by preparing, e.g., three ultracold 
fermionic atoms in an optical potential with four local min¬ 
ima located at the vertices of a square, the results found in our 
work can be checked. This only requires a measurement of 
the 1-particle occupation numbers as performed in [|2f| [ 3j|. 



FIG. 11. Left side: A laser wave (blue) couples to the magnetic 
moment of a system E. However, since E is prepared in a pinned 
state its linear response (right) may be suppressed and would lead to 
a transparency effect (see text). 


Besides such a direct verification of pinning by measuring 
NON an indirect and more fault-tolerant approach is also pos¬ 
sible. It makes use of the fact that pinning potentially restricts 
the dynamics of the corresponding system. This is illustrated 
in Fig. [TO] On the left side, the valence electron cannot decay 
to lower lying energy states since those are already occupied. 
Such a kinematical restriction by Pauli’s exclusion principle 
is generalized by GPC to arbitrary ‘points’ A on the poly tope 
boundary. There, as presented on the right side of Fig. [TOj 
pinned A can never leave the poly tope. A complete mathe¬ 
matical analysis reveals even stronger implications. The lin¬ 
ear response of pinned A to a perturbation of the correspond¬ 
ing 7V-fermion quantum state | ) is restricted to the polytope 
facet & This means that A in leading order cannot leave the 
polytope facet under any small perturbation of the system and 
in particular cannot move to the interior of V. For quasipin¬ 
ning, this holds approximately, as well. In Fig. |TT] a set-up 
is proposed for the indirect experimental verification of pin¬ 
ning by referring to its robustness. A physical system E, e.g., 
the 3-fermion Hubbard model, is perturbed by a laser wave 
(shown in blue) which couples to the magnetic moment M of 
E. By describing the interaction process from the 1-particle 
viewpoint (right side) the linear response of E may signifi¬ 
cantly depend, e.g., on the polarization or the direction of the 
laser wave. In some specific cases (shown in red), the linear 
response of A is expected to be perpendicular to the polytope 
facet. However, this is in contradiction to the existence of 
GPC and therefore the interaction between E and the laser 
wave has to be suppressed. A detector as shown in the set-up 
on the left side of Fig. [TT] can then detect this non-scattered 
laser wave. Furthermore, if E is given by the few-site Hub¬ 
bard model, changing its on-site interaction would allow to 
verify the pinning-nonpinning transitions found in our work. 

By considering not only a single few-site Hubbard model 
but macroscopically many, a GPC-induced transparency- 
effect for that material/gas could be observed. All laser waves 
with a specific polarization, frequency, etc. are not scattered 
and the material is not susceptible to them. 


VII. SUMMARY AND CONCLUSION 


The fermionic exchange statistics does not only imply 
Pauli’s exclusion principle but leads to even stronger restric¬ 
tions on fermionic natural occupation numbers (NON). There¬ 
fore, these so-called generalized Pauli constraints (GPC) have 
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additional influence on the behavior and the properties of 
many-fermion systems. A particular feature is their saturation, 
i.e. the pinning of NON to the boundary of the allowed region. 
We have argued that this pinning effect is very rare, or even 
imp ossible, and that instead quasipinning occurs MH- 

ElESSH. 

As discussed above pinning has strong implications. For 
instance, it leads to a significant simplification of the many- 
fermion quantum state |\h), accompanied by a reduction of 
its entanglement. Due to these strong implications of pinning 
it is highly desirable to find models for which pinning really 
occurs and which can be realized experimentally. 

In form of the few-site Hubbard model we succeeded in 
providing a first system satisfying these requirements. For 
all cases (TV, r) with TV < 5 fermions on r < 5 sites we 
found for many eigenstates finite or even semi-infinite inter¬ 
vals for the relative onsite-interaction u , where their NON 
A (u) are pinned. The end points of those intervals define 
‘critical’ points where sharp pinning-nonpinning or pinning- 
quasipinning transitions are found. As an ‘order parameter’ 
the minimal distance D m i n (u) of A (u) to the polytope bound¬ 
ary dV has been used. We also studied the origin of those 
transitions which turns out to be either a crossing of eigenen- 
ergies or a crossing of NON A i(u), 

The basic reason behind our findings is the presence of ro¬ 
tational (in spin space), and particularly of translational sym¬ 
metry. This immediately implies that the natural spin orbitals 
(the eigenstates of the 1 -particle reduced density operator) are 
the Bloch states | k) multiplied by a spin state | t), | |). Hence 
the natural spin orbitals are known from the very beginning, in 
contrast to atomic and molecular systems. As a consequence 
the NON gain direct physical relevance. 

We have suggested two different experimental set-ups for 
the verification of the pinning-effect, as well as the pinning- 
to-nonpinning transition. Note that the existence of the non- 
pinned ‘phase’ proves that symmetry does not automatically 
implies pinning. Since GPC emerge from the antisymmetry, 
the experimental verification of the pinning effect would also 
provide additional evidence for the antisymmetry of fermionic 
wave functions. 

The focus of the present contribution has been on pinning. 
It will be a challenge for future work to investigate its implica¬ 
tions. One of them is the simplification of the corresponding 
many-fermion quantum state. It will be interesting to study 


how far the simplification of | ^) influences its entanglement 
and the low-lying excitations, responsible for the low temper¬ 
ature behavior. 
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Appendix A: Diagonalization of the subspace Hamiltonian © 

In this section we solve the eigenvalue problems (ITOl) where 
the Hamiltonian is given by Eq. ©. 

The eigenvalues are the roots of the characteristic polyno¬ 
mial P U {E) of ©, i.e. we need to solve the cubic equation 

P U (E) =6u + E(u 2 -9) -2E 2 u + E 3 = 0. (Al) 

Since Appendix[B]requires details on the structure of the roots 
of cubic equations we present in this section a few more de¬ 
tails (see e.g. also |[35|]) than necessary . For the canonic form 

x 3 + ax 2 + for + c = 0 (A2) 

we first consider the quantities Q and R , defined as 

_ a 2 — 3 b _ 2a 3 — 9 ab + 27c 

Q = ,R= - - -. (A3) 

By defining 

0 = arccos (r/VQ*) (A4) 

the three real roots of (1A2I) follow as 

Xk — —2 \[Q cos + g 27rfc ^ - ^ , k = 0,1,2. (A5) 

For our concrete characteristic polynomial (1A11) we have 

'll 2 VL '^ 

Q = - + 3 , R « - . (A6) 

and as real roots (energy eigenvalues) Ei(u),E 2 (u) and 
Es(u) we find 



— \/ u 2 + 27 cos 


1 

- arccos 
3 



27r(l — k ) 
3 


}■ 


Here, we ordered the eigenvalues differently than in (1A51) to (ay (u), /3j {u)^j ( u)),j = 1, 2, 3 of matrix © follow directly 

be consistent with the labeling in Fig. [T] Notice also the sym¬ 
metries of these energies, 

Ei(-u) = -E 3 (u) and E 2 (-u) = -E 2 (u) . (A7) 

The corresponding unnormalized eigenvectors Vj(u) = 
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from ©. We find 

OLj(u) = u + 3 — Ej(u) 
f3j (u) = u — 3 — Ej ( u ) 

7 j(u) = u — 4 Ej(u) + — [ Ej(u ) 2 — 9] . (A8) 


which after division by y | olj (u) | 2 + | j3j (u) | 2 + \^j (u) | 2 

yields the normalized coefficients of the eigenstates (fl 2 l) . 


Appendix B: Calculation of the NON for the superposition ( H9l 

In this appendix we calculate the NON for the more general 
ground state dm Since (IT9l) is not symmetry-adapted to the 
lattice translation anymore, the corresponding 1 -particle re¬ 
duced density operator pi represented w.r.t. the Bloch states 
{| k,a)} is not diagonal anymore. However, since (fl9l) is 
still an 5^-eigenstate, the new natural orbitals are still spin 
eigenstates w.r.t. z-axis and pi splits into blocks according to 
Eq. (IT3l) . The specific normalizations of both blocks together 
with the three equalities CD imply again that for each eigen¬ 
value nj of p\ there exists a corresponding eigenvalue n\ of 
Pi such that their sum equals one, nj + n\_- = 1, j = 1,2,3, 

where the triples {nj} and {nj} are ordered decreasingly. 
Hence, it is sufficient for our purpose to diagonalize just the 
block p\. By representing it w.r.t. the states {\k |)} we obtain 


the following density matrix 

/ \a\ 2 -CC«*7 -CC*a*7 \ 

p\= -eca 7 * IC| 2 |/3| 2 + I?| 2 |7| 2 -rc|/3| 2 

V-cc«7* -eci/3i 2 leiw + icm 2 / 

(Bl) 

The cubic characteristic polynomial p u ( A) reads 

Pu{ A) = + C 2 A 2 + Ci A + Co = 0 (B 2 ) 

c 2 = -1 

Cl = Hu )\ 2 \/ 3 ( u )\ 2 + | a («)| 2 | 7(«)| 2 + \ P ( u )\ 2 I 7 WI 2 
_ (2 — 3 | 7 ( u )| 2 ) | 7 ( m )| 2 | CI 2 1£| 2 
c 0 = -\a(u)\ 2 \f3(u)\ 2 Mu)\ 2 

x (i- 3 |ci 2 i?i 2 -c 3 4 5 6 7 (n 3 -(c*) 3 e 3 ) • 

Its roots follow immediately from Eqs. (1A2K (IA3k (IA4b and 
(IA5I ) in Appendix lAl 

As a consistency check, notice also that according to 
Eq. (IA5b we have in general 

Xi + X 2 + %3 = — C 2 (B3) 

which takes here the value — C 2 = 1 = Tr [p\]. 

Finally, we comment on £, ^-symmetries of p\. First, note 
that the vector of NON fulfills 

A(ix;C ,0 = A(ix;e,C), ( B4 ) 

which follows directly from the invariance of © w.r.t. in¬ 
versions in the Brillouin zone and the form (fl9l) . Moreover, 
blocks pi(u; fulfill 

£,Ce^) = U'p\(u£, QU (B5) 

with the unitary matrix U = diag(l, e z l , e £ l ). The same 
actually also holds for the second block, p[, which is not 
stated here. Thus, 

A(^,Ce^) = A(^,0- (B 6 ) 
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